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======================================================================

*** Introduction

In this section, we will be focussing on conservation of momentum.  We
will be spending a great deal of time studying collisions.

A collision between two objects is defined as an interaction in which
the two objects exert rather large forces on each other over a
relatively short period of time.

Ex: Two erasers colliding in midair [demo 180]
dt on the order of ms -- forces on the order of 100's x weight

More can be said about the initial and final velocities than about
the interaction forces itself.  This makes collisions well-suited
for study using conservation principles, such as conservation of
momentum.

*** Impulse-momentum Theorem
**** Derivation (general) and discussion

F_net = ma
F_net,avg = m a_avg = m Dv/Dt
J_i = F_i,avg Dt; J_net = sum_i J_i
==> J_net = m Dv = Dp  [impulse-momentum theorem]

Alt:
F_net = ma = dp/dt
J_net = int(t_i->t_f) F_net dt = int(t_i->t_f) dp/dt dt = Dp
Note: J_i = int(t_i->t_f) F_i dt


What is impulse?
  F=1000N x-hat, dt=2ms ==> J = Fdt = 2 Ns x-hat
to double impulse...
  F=2000N x-hat, dt=2ms ==> J = Fdt = 4 Ns x-hat
OR
  F=1000N x-hat, dt=4ms ==> J = Fdt = 4 Ns x-hat

Impulse measures effect of force on motion, and depends on force and
how long it acts.


What is momentum?
  m=140g, v=40m/s x-hat (fast-ball) ==> p=mv = 5.6 kg m/s x-hat
to double momentum...
  m=280g, v=40m/s x-hat ==> p = 11.2 kg m/s x-hat
  m=140g, v=80m/s x-hat ==> p = 11.2 kg m/s x-hat

140g baseball starts from rest and is accelerated by force to 40m/s
==> J_net = dp = (5.6 kg m/s x-hat) - (0 kg m/s) = 5.6 Ns x-hat

Momentum equals impulse required to bring object up to speed from rest.

**** Collision examples [140g baseball, direct hit (1D) and foul tip (2D)]

Ball hit by bat (direct); focus on ball  [demo 190a]
  m = 140g = 0.140kg
  vi = -40 m/s x-hat
  vf = +50 m/s x-hat
[before/during/after picture]

Impulse delivered by bat?
  J_net = p_f - p_i = mv_f - mv_i = +12.6 Ns x-hat

Suppose dt=2.0ms.  F_avg?
  F_avg = J/dt = [+12.6 Ns x-hat] / 2.0ms = +6300 N x-hat  [about 1400 lb]
  a_avg = F_avg/m = 4.5x10^3 m/s^2 = 4600(g)

Expected F vs t graph -- 6300N is the average!

Is J_net = J_bat?  Why can we neglect gravity?


Ball hit by bat (glancing); focus on ball  [demo 190b]
  m = 140g = 0.140 kg
  vi = -40 m/s x-hat
  vf = 40 m/s 12deg above -x direction
     = (-39.13 x-hat + 8.32 y-hat) m/s
[before/during/after picture]

Impulse delivered by bat?
  J_net = p_f - p_i  (vectors)

demonstrate -- indicate direction

SI units    x        y
    p_f  -5.478   +1.164
 - [p_i  -5.600    0.000]
-------------------------
  J_net  +0.122   +1.164

J_net = 1.171 Ns, 84 deg above +x direction

Note: 84 deg is exact: J_net bisects angle between v_f and -v_i, due
to equal magnitudes of v_i and v_f.

======================================================================

*** Two-body collisions

**** Conservation of momentum in collisions

Before collision: m_1, m_2, v_1i, v_2i
During collision: F_12, F_21  No external forces!!
After collision: v_1f, v_2f

  J_1,net = J_12 = p_1f - p_1i
+ J_2,net = J_21 = p_2f - p_2i
------------------------------
             0   = p_1f - p_1i + p_2f - p_2i
==> p_1i + p_1f = p_2i + p_2f

Note role of Newton's 3rd law: F_21 = -F_12, so J_21 = -J_12
Also note assumption: no external forces (or neglect them)

Is KE conserved?
  K_1i + K_2i =?= K_1f + K_2f

Not necessarily.  W_12 and W_21 may not cancel because the
displacements are different for the two colliding objects.  More will
said later.

**** 1-D Car examples [stick together / v1f given]

Initial situation (two cars):
  m_1 = 1500 kg
  v_1i = +14 m/s   (driver is texting)
  m_2 = 2000 kg
  v_2i = 0  (stopped at traffic light)


Example 1: stick together after collision
  v_1f = v_2f = v_f

m_1 v_1i + m_2 v_2i = (m_1 + m_2) v_f
+21000 kg m/s = (3500 kg) v_f
v_f = +6.0 m/s

Kinetic energy (kJ):
            1        2      total
initial   147.0      0      147.0
final      27.0    36.0      63.0
change   -120.0   +36.0     -94.0


Example 2: v_1f = +3.0 m/s

m_1 v_1i + m_2 v_2i = m_1 v_1f + m_2 v_2f
+21000 kg m/s = (+4500 kg m/s) + (2000 kg)v_2f
v_2f = +8.25 m/s

Kinetic energy (kJ):
            1        2      total
initial   147.0      0      147.0
final       6.75   68.06     74.81
change   -140.25  +68.06    -72.19


Second collision involves larger momentum transfer.
[demo 210]

KE is lost in both collisions -- less so in second collision.  Other
forms of energy must be involved in such collisions (dominant form
is thermal energy -- usually produced by frictional forces).

======================================================================

**** 2-D example

Initial situation:   [draw picture]
  m_1 = 4.0 kg
  v_1i = 3.0 m/s, due E   [|p_1i| = 12.0 kg m/s]
  m_2 = 3.0 kg
  v_2i = 2.4 m/s, 30 deg W of N  [|p_2i| = 7.2 kg m/s]


Example 3: stick together after collision
  m_1 v_1i + m_2 v_2i = m_1 v_1f + m_2 v_2f = (m_1+m_2) v_f
  v_f = (m_1 v_1i + m_2 v_2i)/(m_1+m_2)

SI units assumed...
              x (E)  y (N)
  m_1 v_1i   +12.00    0
+ m_2 v_2i    -3.60  +6.24
--------------------------
P_tot         +8.40  +6.24    10.46 kg m/s, 36.6 deg N of E
v_f           +1.20  +0.891    1.49 m/s, 36.6 deg N of E

Kinetic energy (J):
            1      2      total
initial   18.00   8.64    26.64
final      4.47   3.35     7.82
change   -13.53  -5.29   -18.82


Example 4: blocks separate and
  v_1f = 1.2 m/s, 57 deg N of E  [|p_1f| = 4.8 kg m/s]
==>
  v_2f = [m_1 v_1i + m_2 v_2i - m_1 v_1f]/m_2

               x (E)   y (N)
   m_1 v_1i   +12.00    0
+  m_2 v_2i    -3.60  +6.24
- [m_1 v_1f    +2.61  +4.03]
----------------------------
m_2 v_2f       +5.79  +2.21    6.19 kg m/s, 20.9 deg N of E
v_2f           +1.93  +0.737   2.06 m/s, 20.9 deg N of E

Kinetic energy (J):
            1      2      total
initial   18.00   8.64    26.64
final      2.88   6.39     9.27
change   -15.12  -2.25   -17.37

======================================================================

**** Launch example

Example 5:
  m_1 = 6800 kg  (cannon on skis)
  m_2 = 34 kg  (ball)
  v_1i = v_2i = 0
  v_2f = +120 m/s
  v_1f = ?  (recoil)

m_1 v_1i + m_2 v_2i = m_1 v_1f + m_2 v_2f
 0 + 0 = m_1 v_1f + m_2 v_2f
==> v_1f = (-m_2/m_1) v_2f = (-34/6800)(+120 m/s) = -0.6 m/s

Kinetic energy (kJ):
            1       2       total
initial    0.0     0.0       0.0
final      1.2   244.8     246.0
change    +1.2  +244.8    +246.0

======================================================================

**** Kinetic energy and collisions

This would be a good time to discuss kinetic energy in each of the
five examples above.

======================================================================

*** Collisions in the center-of-mass reference frame
**** Collisions viewed in the center of mass frame

[See "Two object systems" in sys.txt]

vcm_i = [m1 v1i + m2 v2i]/M
vcm_f = [m1 v1f + m2 v2f]/M
where M = m1 + m2

Conservation of momentum ==> vcm_i = vcm_f.  This makes sense, since
F_ext = 0 for the two-particle system during the collision.  In fact,
this is usually how momentum conservation is derived.

The center-of-mass reference frame is defined by placing the origin
at the center of mass for all time.  This implies that vcm = 0 for
all time.

Since vcm=0 in the center-of-mass reference frame,
  v1i = +(m2/M) v12i ; v1f = +(m2/M) v12f
  v2i = -(m1/M) v12i ; v2f = -(m1/M) v12f

[draw picture of collision in this frame, and show v_12i and v_12f]

This implies that the collision can be completely characterized by
v_12i and v_12f.  Any values of v_12i and v_12f would be consistent
with momentum conservation.

It should be noted that ANY collision can be viewed in its center
of mass reference frame.


Car examples (1 and 2):

m1 = 1500 kg, v1i = +14 m/s
m2 = 2000 kg; v2i = 0

M = 3500 kg
vcm = [m1 v1i + m2 v2i]/M = +6.0 m/s
v12i = v1i - v2i = +14 m/s

Note that vcm is the final velocity that would result if the
colliding objects stick together, even if that is not what really
happens.

Converting each velocity to the cm frame yields:
  v1ci = v1i - vcm = +8.0 m/s
  v2ci = v2i - vcm = -6.0 m/s    [draw picture]

For example 1, the cars stick together so v1f = v2f = vf = vcm.  Thus,
v12f = 0 and v1cf = v2cf = 0.  There is no relative motion after the
collision.

For example 2, v1f was given (+3.0m/s) and v2f was calculated from
momentum conservation (+8.25m/s), and so
  v12f = v1f - v2f = -5.25m/s
Converting each final velocity to the cm frame yields:
  v1cf = v1f - vcm = -3.0 m/s
  v2cf = v2f - vcm = +2.25 m/s

Note that v_12i and v_12f can be calculated directly without first
converting to the cm reference frame.


2-D examples (3 and 4):

Initial situation:
  m_1 = 4.0 kg
  v_1i = 3.0 m/s, due E
  m_2 = 3.0 kg
  v_2i = 2.4 m/s, 30 deg W of N
Example 3: stick together after collision
Example 4: blocks separate and v_1f = 1.2 m/s, 57 deg N of E

All velocities in m/s...
            x(E)     y(N)     magitude  direction (ccw from E)
Initial...
v1i        +3.000    0.000      3.000       0.0
v2i        -1.200   +2.078      2.400     120.0
vcm        +1.200   +0.891      1.494      36.6
v12i       +4.200   -2.078      4.686     333.7
v1ci       +1.800   -0.891      2.008     333.7
v2ci       -2.400   +1.188      2.678     153.7

Example 3...
v1f=v2f    +1.200   +0.891      1.494      36.6
v12f        0.000    0.000      0.000     -----
v1cf=v2cf   0.000    0.000      0.000     -----

Example 4...
v1f        +0.654   +1.006      1.200      57.0
v2f        +1.929   +0.737      2.064      20.9
v12f       -1.275   +0.270      1.303     168.1
v1cf       -0.546   +0.116      0.559     168.1
v2cf       +0.729   -0.154      0.745     348.1

Alternative method of calculating v2f in example 4:
(1) v1cf = v1f - vcm
(2) v12f = (+M/m2)v1cf
(3) v2cf = (-m1/M)v12f = (-m1/m2)v1cf
(4) v2f = v2cf + vcm

Picture arrangements:
            lab frame      c.m. frame
initial
final (1)
final (2)

For initial/lab, indicate v1i, v2i, and vcm
For initial/cm, indicate v1ci, v2ci, and v12i
For final(1)/lab, indicate vf=v1f=v2f
For final(1)/cm, both objects are at rest
For final(2)/lab, indicate v1f and v2f
For final(2)/cm, indicate v1cf, v2cf, and v12cf


Example 5 (cannon launch):
  m_1 = 6800 kg  (cannon on skis)
  m_2 = 34 kg  (ball)
  v_1i = v_2i = 0
  v_2f = +120 m/s
Momentum conservation led to recoil velocity...
  v_1f = -0.6 m/s

Since v_1i and v_2i are both zero, v_cm = 0, and so this "collision"
is already being viewed in the center-of-mass frame.  There is nothing
to "convert".

Note that v_12i=0 (no relative motion to start), but that
v_12f = v_1f - v_2f = -120.6m/s as a result of the explosive launch.

======================================================================

**** The elastic parameter

Given that v_12i and v_12f completely determine the collision, it
follows that collisions can be completely characterized by a
comparison between v_12i and v_12f in terms of magnitude and
direction.

The elastic parameter (for lack of a better term) is determined by
comparing magnitudes...
  |v_12f| = alpha |v_12i|

alpha is a measure of the "liveliness" of a collision.

It can also be related to gain/loss in kinetic energy during the
collision.

alpha = 0 (|v_12f|=0)  ==> completely inelastic (stick together)
alpha < 1 (|v_12f|<|v_12i|) ==> (partially) inelastic -- frictional losses (MOST COMMON)
alpha = 1 (|v_12f|=|v_12i|) ==> elastic collision (spring-like, conservative)
alpha > 1 (|v_12f|>|v_12i|) ==> super-elastic / explosive (energy source required)

There is a related parameter called the "coefficient of restitution",
which we will introduce later.


Kinetic energies:

In the center-of-mass frame...

v1 = (+m2/M)v12
v2 = (-m1/M)v12

K_sys = (1/2)m1 v1^2 + (1/2)m2 v2^2
  = (1/2)m1 (m2/M v12)^2 + (1/2)m2 (m1/M v12)^2
  = (1/2)[m1 m2^2/M^2 + m1^2 m2/M^2] v12^2
  = (1/2)(m1m2/M^2)(m1+m2) v12^2
  = (1/2)(m1m2/M) v12^2
  = (1/2)mu v12^2
where mu = (1/m1 + 1/m2)^-1 = m1m2/M is the reduced mass of the
system.

Thus,
K_sys_f = (1/2)mu v_12f^2
  = (1/2)mu (alpha v_12i)^2
  = alpha^2 (1/2)mu v_12i^2
  = alpha^2 K_sys_i

In a general reference frame...

v1 = vcm + (m2/M)v12
v2 = vcm - (m1/M)v12

K_sys = (1/2)m1 |vcm + (m2/M)v12|^2 + (1/2)m2 |vcm - (m1/M)v12|^2
  =   (1/2)m1[vcm^2 + (m2/M)^2 v12^2 + 2(m2/M)vcm dot v12]
    + (1/2)m2[vcm^2 + (m1/M)^2 v12^2 - 2(m1/M)vcm dot v12]
  = (1/2)[m1+m2]vcm^2 + (1/2)[m1m2^2/M^2 + m2m1^2/M^2]v12^2
    + (1/2)[2m1m2/M - 2m2m1/M](vcm dot v12)
  = (1/2)M vcm^2 + (1/2)mu v12^2   (third term zero)
    ^^^ K_cm ^^^   ^^^ K_rel ^^^

During a collision, v_cm never changes, and so K_cm does not
change.  However, K_rel (and thus K_sys) can change due to changes
in relative motion (|v12|).  In fact, K_rel will change by a factor
of alpha^2.

This implies the following:
  alpha=1 (elastic) ==> KE is conserved
  alpha<1 (inelastic) ==> KE decreases
  alpha>1 (super-elastic) ==> KE increases

======================================================================

**** The deflection parameter

The deflection parameter (for lack of a better term) is determined by
comparing the directions of v_12f and v_12i.

The deflection parameter is defined as the angle, phi, between v_12f
and v_12i.  phi is also the angle between v1cf and v1ci, as well as
v2cf and v2ci.   [draw picture showing collision in cm frame]

This angle can vary from phi=0deg to phi=180deg, and is related to how
"direct" or "head-on" the collision is.  It can also be related to the
impulse delivered by the collision forces during the collision.

Small angles indicate very little deflection, and represent collisions
with minimal impulse.

Note: phi=0 represents either a complete miss (alpha=1) or a
"pass-through collision", where the objects slow down (if alpha<1) as
they pass through each other (e.g. bullet through a block).

Larger angles represent more direct collisions and give rise to larger
impulses.

phi=180deg corresponds to a direct head-on collision, and provides for
maximum impulse.

There is a related parameter called the "impact parameter", which is
defined in terms of a "distance of closest approach" of the two
objects, if their interactions were turned off [draw picture].  Larger
impact parameters tend to lead to smaller deflection angles as a
general rule.


Impulses:

J12 = dp1 = m1(v1f - v1i) = m1(+m2/M)(v12f - v12i) = mu(v12f - v12i)
J21 = dp2 = m2(v2f - v2i) = m2(-m1/M)(v12f - v12i) = -mu(v12f - v12i)

[draw the subtraction triangle v12f-v12i]

Thus,
|J|^2 = mu^2 |v12f - v12i|^2
  = mu^2 (v12f - v12i) dot (v12f - v12i)
  = mu^2 v12f^2 - 2 v12f dot v12i + v12i^2
  = mu^2 v12i^2 [alpha^2 - 2 alpha cos(phi) + 1]

It follows that
  |J| = mu |v12i| sqrt[alpha^2 - 2 alpha cos(phi) + 1]
For fixed alpha, |J| increases as phi increases from 0 deg to 180 deg.

It is more interesting if we fix phi and allow alpha to vary:

alpha^2 - 2 alpha cos(phi) + 1
  = alpha^2 - 2 alpha cos(phi) + cos(phi)^2 + sin(phi)^2
  = (alpha - cos(phi))^2 + sin(phi)^2
and so
  |J| = mu |v12i| sqrt[(alpha - cos(phi))^2 + sin(phi)^2]

If phi >= 90 deg, |J| increases with increasing alpha.

If phi < 90 deg, |J| will decrease as alpha is increased from alpha=0
to alpha=cos(phi) (minimum J), and then |J| will increase thereafter.

======================================================================

**** Examples revisited

Example 1:
  v12i = +14 m/s
  v12f = 0
==> alpha = 0, phi undefined  (completely inelastic)
==> e = 0, phi' = 0 deg, phi'' undefined

Example 2:
  v12i = +14 m/s
  v12f = -5.25 m/s
==> alpha = 0.375 (3/8), phi = 180 deg  (partially inelastic; head-on)
==> e = 0.375, phi' = phi'' = 0 deg

Example 3:
  v12i = 4.686 m/s, 333.7 deg ccw +x
  v12f = 0
==> alpha = 0, phi undefined  (completely inelastic)
==> e = 0, phi' = 0 deg, phi'' undefined

Example 4:
  v12i = 4.686 m/s, 333.7 deg ccw +x
  v12f = 1.303 m/s, 168.1 deg ccw +x
==> alpha = 0.278, phi = 165.6 deg  (partially inelastic, nearly head-on)
==> e = 0.273, phi' = 3.12 deg, phi'' = 11.29 deg

Note that none of these collisions are elastic, even though two of
them involve the colliding objects separating after the collision.
Some KE is lost in all cases.

It is difficult to construct an elastic collision "by accident".

Example 5: cannon launch
  v12i = 0
  v12f = -120.6 m/s
==> alpha = infinity, phi undefined  (explosive -- literally)
==> e = infinity, phi' undefined, phi'' = 0 deg
======================================================================

**** The angle of incidence and coefficient of restitution

We have already learned about describing collisions using (alpha,phi).
Now we will define another equivalent set of parameters that can be
used to characterize collisions, namely (e,phi'), where
  e = coefficient of restitution
  phi' = angle of incidence

It will turn out that (alpha,phi) are more convenient to use in
calculations, since these parameters relate v12i and v12f directly.
This will allow us to quickly calculate final velocities from initial
velocities, etc.

However, (e,phi') may be more physically relevant, since they are
defined in terms of the direction of the impulse.  In any case, we
will discover later how to convert between them.

[draw picture showing colliding objects in cm view:
(1) objects are in contact
(2) initial and final velocities are shown for each (labeled for 1)
(3) dotted line used for plane of collision
(4) n-hat represents direction of interaction force (and thus J12)]

The impulse is given by
  J12 = mu(v12f - v12i)
Its direction is given by the unit vector
  n-hat = J12 / |J12|

For a surface contact collision without significant friction, we
expect n-hat to be perpendicular to the surface at the point of
contact between the colliding objects.  The collision forces will
largely be normal (compression) forces involving both objects.

[now draw v12i, v12f, dv12 triangle -- indicate phi]

The angle of incidence, phi', is defined to be the angle between
-v12i and n-hat (J12).  [indicate in diagram -- difference between
phi and phi' should be evident]

We can also define the angle of reflection, phi'', to be the angle
between J12 and +v12f.  [indicate in diagram]

If we imagine a collision between a ball and a "large" surface,
such as a wall/floor (attached to a planet), the collision can
be viewed more simply...

[now draw ball bouncing off of floor; indicate v1i, v1f, and n-hat
straight upward and indicate phi' and phi'']

If we split the relative velocities into components parallel to
and perpendicular to n-hat:
  v12i = v12i_par + v12i_perp   [indicate this in ball-floor diagram]
  v12f = v12f_par + v12f_perp
we find that v12f_perp = v12i_perp -- v12_perp does not change.
That is because there is no impulse in that direction.

However, v12_par DOES change -- in fact, it often reverses direction.
If v12i_par = -v12f_par exactly, then |v12i|=|v12f| overall, and the
collision would be elastic.  However, that doesn't always happen.

So we define the coefficient of restitution to represent the
actual relationship:

e = v12f_par / (-v12i_par)
  = (v12f dot n-hat) / (-v12i dot n-hat)

We expect this coefficient to represent, in some sense, the efficiency
of the normal force to restore the speed of the object in the
direction of the normal force.  As such, it is not expected to depend
(much) on phi', although it will certainly depend on the physical
properties of the colliding objects themselves (contrast with alpha).

* e = 1 would represent an elastic collision.
* e > 1 would represent a (normal) super-elastic collision.
* 0<=e<1 is expected for most collisions with some KE loss.
* e < 0 can happen for certain low deflection collisions which are
significantly lossy (in terms of KE).
* e < -1 is extremely rare, and generally requires an impulse that
actually enhances v12_par in the direction it is already pointing in
(an explosion "from behind", as the objects pass through each other).

What if there is significant friction between the two colliding
objects (resulting in a component of impulse parallel to the contact
surface)?  Or what if the collision itself is not a surface contact
collision at all (e.g., two protons fired towards each other and
scattering due to Coulomb repulsion)?

In that case, it may make sense to continue defining n-hat in terms of
the direction of J12 (never mind the nature of the forces involved),
in which case the relationships we develop between (e,phi') and
(alpha,phi) (for example) would remain unchanged.

However, e and phi' may lose their physical meaning -- e would no
longer be based strictly on compression/restoration efficiency, and
thus e may depend strongly on phi', and not just on the objects
involved in the collision.

Perhaps, to obtain better modelling, it may be more prudent to
continue to define n-hat to be the contact surface normal and
take friction into account explicitly, perhaps with additional
parameters that define how v12_perp changes during the collision.

I was unable to resolve this with research -- most online resources
which discuss the coefficient of restitution ignore this issue
altogether -- it is assumed that J12 was perpendicular to the contact
surface.

======================================================================

**** Conversions

[Draw ball bouncing off floor picture, indicate v1i, v1f, and n-hat
straight upward and indicate phi, phi', and phi'']

[Now draw v12i-v12f-dv12 triangle with dv12 straight up.  Draw
perpendicular (horizontal) from v12i-v12f vertex to dv12.  Indicate
angles (phi, phi', phi'') and components v12_perp, v12i_par,
v12f_par.]

The impulse J12 = mu (v12f - v12i), and is in the direction of
dv12.  n-hat will a unit vector in that direction.

(e,phi') -> (alpha,phi):

v12f_perp = v12i_perp = v12i sin(phi')
v12f_par = -e v12i_par = +e v12i cos(phi')
|v12f| = v12i sqrt[sin(phi')^2 + e^2 cos(phi')^2]
==> alpha = sqrt[sin(phi')^2 + e^2 cos(phi')^2]
phi'' = arctan(v12f_perp / v12f_par)
  = arctan( [v12i sin(phi')] / [e v12i cos(phi')] )
  = arctan((1/e)tan(phi'))
==> phi = 180deg - phi' - phi''
        = 180deg - phi' - arctan((1/e)tan(phi'))

(alpha,phi) -> (e,phi'):

v12i dot v12i = v12i^2 = v12i^2 [1]
v12i dot v12f = |v12i| |v12f| cos(phi) = v12i^2 [alpha cos(phi)]
v12f dot v12f = v12f^2 = v12i^2 [alpha^2]

e = [v12f dot n-hat] / [-v12i dot n-hat]
  = [v12f dot J12] / [-v12i dot J12]
  = mu v12i^2 [alpha^2 - alpha cos(phi)] / mu v12i^2 [1 - alpha cos(phi)
  = [alpha^2 - alpha cos(phi)] / [1 - alpha cos(phi)]
  = alpha [alpha - cos(phi)] / [1 - alpha cos(phi)]

-v12i dot J12 = mu v12i^2 [1 - alpha cos(phi)]
  = |v12i| |J12| cos(phi')

Recall |J12| = mu v12i sqrt[alpha^2 - 2 alpha cos(phi) + 1]

cos(phi') = mu v12i^2 [1 - alpha cos(phi)] / mu v12i^2 sqrt[alpha^2 - 2 alpha cos(phi) + 1]
  = [1 - alpha cos(phi)] / sqrt[alpha^2 - 2 alpha cos(phi) + 1]

Summary:

  phi'' = arctan[(1/e)tan(phi')]

  alpha = sqrt(sin(phi')^2 + e^2 cos(phi')^2]
  phi = 180deg - phi' - phi''

  e = alpha [alpha - cos(phi)] / [1 - alpha cos(phi)]
  cos(phi') = [1 - alpha cos(phi)] / sqrt[alpha^2 - 2 alpha cos(phi) + 1]

Some things to note:

* alpha = 1 if and only if e = 1.
* alpha = 0 if and only e = 0 and phi' = 0deg
* either |e| <= alpha <= 1 or |e| >= alpha >= 1
* phi = 180 deg ==> phi' = 0deg and e = alpha
* phi = 0deg ==> phi' = 0 deg (alpha<1) or 180 deg (alpha>1) and e = -alpha
* if alpha = 1, then n-hat bisects the angle between -v12i and v12f:
    phi'' = phi'
    phi = 180deg - 2 phi'
    phi' = (1/2)(180deg - phi)
 * for e>0, as phi decreases from 180 deg to 0 deg, phi' will increase
from 0 deg to 90 deg, and alpha will increase from e to 1.
* if phi > 90 deg, then e > 0
* if 1 >= alpha > cos(phi), then e > 0
* if alpha < cos(phi), then e < 0

Revisiting example 4:
  alpha = 0.278
  phi = 165.6 deg
Thus
  e = 0.273
  phi' = 3.12 deg

======================================================================

**** Measuring the coefficient of restitution (demo)

[projectile diagram]

For a ball either dropped (v0x=0) or launched horizontally (v0x>0)
from a height h0 above the ground, the velocity of the object just
before it hits the ground is given by
  v1 = v0x x-hat - sqrt(2*g*h0) y-hat   [projectile equations]

The ball then collides with the ground, after which
  v2x = v1x     [assume minimal friction -- assured by "rolling" launch]
  v2y = -e v1y  [definition of coefficient of restitution]
      = +e sqrt(2*g*h0)
The ball then reaches a maximum height h3 afterwards, where
  h3 = v2y^2/(2*g)
     = e^2 2*g*h0/(2*g)
     = e^2 h0

The coefficient of restitution can be measured by measuring h0 and h3
and calculating
  e = sqrt(h3/h0)

Note that the angle of incidence for the collision will be
  phi' = arctan(v1x/|v1y|)
The time required for the object to fall can be calculated
  dt01 = sqrt(2h0/g)
which allows us to express the initial velocity in terms of the
horizontal distance from launch to the collision:
  v1x = dx01 / sqrt(2h0/g)
It follows that
  v1x/|v1y| = dx01 / [sqrt(2h0/g) sqrt(2*h0*g)]
    = dx01 / [2 h0]
    = [dx01/2] / h0
and so the angle of incidence will be
  phi' = arctan([dx01/2]/h0)

Video demonstration:  [see 215.m4v]

Position measurements on screen:
* Initial release (vertical): y = 12.7cm  [note: +y is down -- thanks ruler]
* Hit floor: y = 27.2cm
* Back to rest: y = 19.3cm
* Vertical contact point: x = 6.0cm
* Contact point (trial 2): x = 10.5cm
* Contact point (trial 3): x = 16.6cm
* Contact point (trial 4): x = 19.4cm

Screen to real-life scale conversion unknown, but always cancels out:

h0 = 27.2cm - 12.7cm = 14.5cm
h3 = 27.2cm - 19.3cm =  7.9cm
e = sqrt(h3/h0) = 0.738 (0.74)
Trial 1:
  phi' = 0deg
  alpha = 0.738
Trial 2:
  dx01 = 10.5cm - 6.0cm = 4.5cm
  phi' = arctan([4.5cm/2]/14.5cm) = 8.8 deg
  alpha = 0.745
Trial 3:
  dx01 = 16.6cm - 6.0cm = 10.6cm
  phi' = arctan([10.6cm/2]/14.5cm) = 20.1 deg
  alpha = 0.774
Trial 4:
  dx01 = 19.4cm - 6.0cm = 13.4cm
  phi' = arctan([13.4cm/2]/14.5cm) = 24.8 deg
  alpha = 0.790

Use horizontal ruler to demonstrate that h3 (and thus e) does not
appear to depend on trial (phi' value).


If alpha had remained constant (0.738) instead of e, then e
would have changed to
  e = sqrt[alpha^2 - sin(phi')^2]/[cos(phi')^2]
    = 0.669
for trial 4, which leads to a height h3 = e^2 h0 = 6.5cm, which is
1.4cm below the h3 value of the first trial.  That would have been
noticeable.

[vertical distance from time indicator bar to top of timer strip on
Quicktime player is about 1.5cm]

======================================================================

*** Collisions involving an initially stationary target
**** General solution

For a general (2-D) collision, if we are given the masses and initial
velocities, momentum conservation will provide us with two equations
for the four unknown quantities (v1fx, v1fy, v2fx, v2fy).

This implies that two additional quantities are needed in order
to solve for the unknown final velocities.

One option is to stipulate that the two objects stick together
after the collision (v1f = v2f).  This provides two additional
equations, which allow us to solve for the common final velocity.

Another option is to also give one of the two final velocities.
This allows us to use momentum conservation to solve for the other
one.

We have done numerous examples of each.

A third option, which we presently consider, is to specify alpha
and phi.  (Alternatively, we could specify e and phi', but in that
case, it is best to convert to alpha and phi.)

We will make the symplifying assumption that one of the two
objects is initially stationary (v2i=0) -- that object is known as
the "target".  We will also impose coordinate choices: the initial
velocity of the "incident" is in the +x direction, and v12f will
be directed phi counter-clockwise from +x in the x-y plane.

Given: v1i (+x), v2i=0, m1, m2, alpha, phi
Calculate: v1f, v2f

Define beta = m2/m1.

v12i = v1i - v2i = v1i = +v1i x-hat
v12f = |v12f| (cos(phi)x-hat + sin(phi)y-hat)
     = alpha v1i (cos(phi)x-hat + sin(phi)y-hat)
vcm = [m1 v1i + m2 v2i]/M
    = (m1/M)v1i
    = v1i/(1 + beta)
    = [|v1i|/(1+beta)] x-hat

Momentum conservation is equivalent to the assertion that vcmf = vcmi,
and so

v1f = vcm + (m2/M)v12f
    = [v1i/(1+beta)] x-hat + [beta/(1+beta)] alpha v1i [cos(phi)x-hat + sin(phi)y-hat]
    = [v1i/(1+beta)] [(1 + alpha beta cos(phi)) x-hat + (alpha beta sin(phi)) y-hat]
      ^^^^^^^^^^^^^^--- |vcm|

v2f = vcm - (m1/M)v12f
    = [v1i/(1+beta)] x-hat - [1/(1+beta)] alpha v1i [cos(phi)x-hat + sin(phi)y-hat]
    = [v1i/(1+beta)] [(1 - alpha cos(phi)) x-hat - (alpha sin(phi)) y-hat]
      ^^^^^^^^^^^^^^--- |vcm|

It is a simple matter to check that momentum is conserved:
  v1f + beta v2f = v1i

[Draw diagram showing vcm, v1fc, v2fc, v1f, v2f
Label the vectors and also indicate the angle phi.]

======================================================================

**** 1-D collisions

Head-on collisions:

phi = 180 deg
cos(phi) = -1, sin(phi) = 0

v1f = [v1i/(1+beta)] [(1 - alpha beta) x-hat]
    = [(1 - alpha beta)/(1 + beta)] v1i
v2f = [v1i/(1+beta)] [(1 + alpha) x-hat]
    = [(1 + alpha)/(1 + beta)] v1i

Confirmation:
v12f = v1f - v2f
     = [(1 - alpha beta - 1 - alpha)/(1 + beta)] v1i
     = [-alpha (1 + beta) / (1 + beta)] v1i
     = (-alpha)v1i
     = (-alpha)v12i

Note that v1f and v2f are in the +/- x direction.

v2f is always in the +x direction (except as beta->infinity, where v2f
is (nearly) zero).

v1f is either +x, -x, or zero, depending on alpha*beta:
  alpha beta > 1  ==> v1f in -x direction (backwards)
  alpha beta = 1  ==> v1f = 0
  alpha beta < 1  ==> v1f in +x direction (forwards)

[diagram showing before and after, with all three possibilities]]

This makes sense: more lively head-on collisions (larger alpha) are
more likely to send the incident backwards, as are collisions
involving a relatively more massive target.

Setting alpha = 1 gives the "famous" equations for 1-D (head-on)
elastic collisions:
  v1f = [(1-beta)/(1+beta)] v1i = [(m1-m2)/(m1+m2)] v1i
  v2f = [(1+1)/(1+beta)] v1i = [2m1/(m1+m2)] v1i

If we consider the limit beta->infinity (e.g., object makes head-on
collision with a wall, connected to earth), we get
  v1f = [(1/beta - alpha)/(1/beta + 1)] v1i (dividing num and den by beta)
      = [(-alpha)/(1)] v1i = (-alpha)v1i
Object will always bounce backwards, with a speed of alpha times the
initial speed.  If the collision is elastic (alpha=1), the object will
bounce off the wall with the same speed it came in.

[diagram showing object bouncing off of a wall]


Pass-through collisions:

phi = 0 deg
cos(phi) = +1, sin(phi) = 0

An obvious example of this type of collision is a bullet being shot
through a block, and passing out the other side.

v1f = [v1i/(1+beta)] [(1 + alpha beta) x-hat]
    = [(1 + alpha beta)/(1 + beta)] v1i
v2f = [v1i/(1+beta)] [(1 - alpha) x-hat]
    = [(1 - alpha)/(1 + beta)] v1i

Confirmation:
v12f = v1f - v2f
     = [(1 + alpha beta - 1 + alpha)/(1 + beta)] v1i
     = [alpha (1 + beta) / (1 + beta)] v1i
     = alpha v1i
     = alpha v12i

Again, v1f and v2f are in the +/- x direction.

v1f is always in the +x direction.

v2f is either +x, -x, or zero, depending on alpha:
  alpha > 1 (super-elastic)  ==> v2f in -x direction (backwards)
  alpha = 1 (elastic/miss)   ==> v2f = 0
  alpha < 1 (inelastic)      ==> v2f in +x direction (forwards)
v2f will also (nearly) equal zero in the limit beta->infinity.

[diagram showing before and after, with all three possibilities]

This makes sense: objects pass through each other, so it makes
sense that object 1 continues forward.  If the collision is
explosive, the explosion will send the target backwards as the
incident passes through it.  An "elastic pass-through" collision
is a complete miss (v1f=v1i, v2f=v2i).  If the collision is
inelastic, the target will get dragged along as the incident
passes through it (presumably the relative velocity between the
two is decreased by friction).

======================================================================

**** 1-D examples

Head-on collisions: Car examples

v1i = (+14 m/s) x-hat
m1 = 1500 kg; m2 = 2000 kg  (beta=4/3)
alpha = 0, phi not important  (completely inelastic)

vcm = v1i/(1+beta) = (+6.0m/s) x-hat
v1f = v2f = vcm = (+6.0m/s) x-hat

-----

same v1i, v2i, m1, m2
alpha = 0.20, phi = 180 deg  (head-on inelastic)

v1f = vcm (1 - alpha beta) x-hat = (+4.4 m/s) x-hat
v2f = vcm (1 + alpha) x-hat = (+7.2 m/s) x-hat

-----

same v1i, v2i, m1, m2
alpha = 1, phi = 180 deg  (head-on elastic)
v1f = vcm (1 - alpha beta) x-hat = (-2.0 m/s) x-hat
v2f = vcm (1 + alpha) x-hat = (+12.0 m/s) x-hat

Note that v12f = v1f - v2f = (-14 m/s) x-hat = -v12i.

You should check momentum conservation in each case, and
compare kinetic energy values in each case.  KE will be
conserved in the last case.

-----

Pass-through collisions: Bullet passing through a wooden block

v1i = (800 m/s) x-hat
m1 = 10.0 g = 0.010 kg, m2 = 10.0 kg  (beta=1000)
alpha = 0, phi not important  (bullet is stuck in block)

vcm = v1i/(1+beta) = (+0.7992 m/s) x-hat
v1f = v2f = vcm = (+0.7992 m/s) x-hat

Virtually all of KE is lost.

-----

same v1i, v2i, m1, m2
alpha = 0.80, phi = 0 deg  (bullet passes through, inelastic)

v1f = vcm (1 + alpha beta) x-hat = (+640.16 m/s) x-hat
v2f = vcm (1 - alpha) x-hat = (+0.16 m/s) x-hat

Some KE lost.

-----

same v1i, v2i, m1, m2
alpha = 1.2, phi = 0 deg  (explosive)

v1f = vcm (1 + alpha beta) x-hat = v1i x-hat = (+959.84 m/s) x-hat
v2f = vcm (1 - alpha) x-hat = (-0.16 m/s) x-hat

KE is actually gained.  Somebody must have packed the block with gun
powder.

======================================================================

**** 2-D examples

Show spreadsheet...

v1i = (4.686 m/s) x-hat
m1 = 4.0kg, m2 = 3.0kg  (beta = 0.75)
alpha = 0.278, phi = 165.6 deg  (inelastic, nearly head-on)

This is example 4, except transformed to a reference frame where m2 is
initially at rest, and the coordinate system is rotated so that v1i is
in the +x direction.  The alpha and phi for this example were
calculated earlier.

vcm = v1i/(1+beta) = (+2.678 m/s) x-hat
v1f = vcm [(1 + alpha beta cos(phi)) x-hat + (alpha beta sin(phi)) y-hat]
    = (2.137 x-hat + 0.139 y-hat) m/s
    = 2.141 m/s, 3.72 deg ccw from +x
v2f = vcm [(1 - alpha cos(phi)) x-hat - (alpha sin(phi)) y-hat]
    = (3.399 x-hat - 0.185 y-hat) m/s
    = 3.404 m/s, -3.12 deg ccw from +x (3.12 deg cw)
                                        ^^^^--- Note: this equals phi'

Note the angle between the final velocities is quite small -- this collision
is nearly head-on, and alpha is quite small.

-----

Same v1i, m1, m2
alpha = 1, phi = 165.6 deg  (elastic version, same phi)

v1f = vcm [(1 + alpha beta cos(phi)) x-hat + (alpha beta sin(phi)) y-hat]
    = (0.733 x-hat + 0.499 y-hat) m/s
    = 0.887 m/s, 34.29 deg ccw from +x
v2f = vcm [(1 - alpha cos(phi)) x-hat - (alpha sin(phi)) y-hat]
    = (5.271 x-hat - 0.666 y-hat) m/s
    = 5.313 m/s, -7.2 deg ccw from +x (7.2 deg clockwise)

The angle between v1f and v2f is now much larger (consistent with increase
in alpha).

Also note that |v2f| is larger, indicating a larger impulse.  For
phi > 90 deg, higher alpha implies greater impulse.

Kinetic energy (J):
            1        2      total
initial   43.92      0      43.92
final      1.57    42.35    43.92
change   -42.35   +42.35      0

-----

Same v1i, m1, m2
alpha = 0.278, phi = 80.0 deg  (original alpha, more glancing)

v1f = vcm [(1 + alpha beta cos(phi)) x-hat + (alpha beta sin(phi)) y-hat]
    = (2.775 x-hat + 0.550 y-hat) m/s
    = 2.829 m/s, 11.21 deg ccw from +x
v2f = vcm [(1 - alpha cos(phi)) x-hat - (alpha sin(phi)) y-hat]
    = (2.548 x-hat - 0.733 y-hat) m/s
    = 2.652 m/s, -16.05 deg ccw from +x (16.05 deg clockwise)

======================================================================

**** Notable features of 2-D collisions

v1f = vcm + (m2/M)v12f
    = [v1i/(1+beta)] [(1 + alpha beta cos(phi)) x-hat + (alpha beta sin(phi)) y-hat]
v2f = vcm - (m1/M)v12f 
    = [v1i/(1+beta)] [(1 - alpha cos(phi)) x-hat - (alpha sin(phi)) y-hat]
beta = m2/m1

[Re-draw diagram showing vcm, v1fc, v2fc, v1f, v2f.
(1) Label the vectors and also indicate the angle phi.
(2) Note the magnitudes of v1fc (alpha*beta*vcm) and v2fc (alpha*vcm).
(3) Also indicate phi1 and phi2 (angles between v1i and v1f/v2f in the
lab frame).]
phi12 = phi1 + phi2 = angle between v1f and v2f in the lab frame

Note that all lengths are proportional to vcm, and so we get the same
angles from the following structure
[new diagram with vcm divided out -- same angles indicated]

Some things to note:

(1) If phi and beta are fixed and alpha is increased, both phi1
and phi2 increase.  This also implies phi12 (angle between v1f
and v2f) increases.  Collisions with increased elastic parameter
"spread out more".

(2) If phi and alpha are fixed and beta is increased, then phi1 is
increased, phi2 is unchanged, and again phi12 increases.  Collisions
with relatively higher mass targets also "spread out more".

(3) If alpha=beta=1 (elastic collision with equal masses) then both
triangles are isosceles:  [redraw them]

    phi1 + phi1 + (180deg-phi) = 180 deg
 +  phi2 + phi2 + phi = 180 deg
-----------------------------------
    2(phi1 + phi2) + 180 deg = 360 deg
==> phi12 = phi1 + phi2 = 90 deg  -- colliding objects leave at 90 deg
angle.  This is a well-known result (pool players use this result all
the time).

(4) It follows that if alpha <= 1 and beta <= 1, then phi12 <= 90 deg;
and also if alpha >= 1 and beta >= 1, then phi12 >= 90 deg.  All bets
are off if alpha < 1 and beta > 1 or vice versa.

We can also fix alpha and beta and consider variations in phi...


Focus on phi1...

[Draw diagram focussing on v1f:
(1) Draw vcm.
(2) Draw semi-circle of radius alpha beta vcm (alpha beta < 1).
(3) Label far point as 0 deg and near point as 180 deg
(4) Draw several vectors v1f -- include phi_1,max]

As phi varies from 0 deg to 180 deg, v1cf spins around a semi-circle
centered at the head of vcm.  v1f will point from the tail of vcm to
the appropriate point on that semi-circle.

alpha beta < 1:

phi: 0 deg -> 180 deg
phi1: 0 deg -> phi_1,max -> 0 deg
where phi_1,max = arcsin(alpha beta)

alpha beta = 1:  [redraw diagram]

phi: 0 deg -> 180 deg
phi1: 0 deg -> 90 deg

Note that v1f->0 as phi->180 deg and phi1->90 deg in this case.

alpha beta > 1:  [redraw diagram]

phi: 0 deg -> 180 deg
phi1: 0 deg -> 180 deg  (all angles, including backward scattering, are possible)

Note that larger scattering angles for object 1 occur as either beta
or alpha increases.  Features of the collision that work in favor of
these larger scattering angles are
(1) More "lively" collision (increased alpha)
(2) More massive target


Focus on phi2...

[Draw diagram focussing on v2f:
(1) Draw vcm.
(2) Draw semi-circle of radius alpha vcm (alpha < 1).
(3) Label near point as 0 deg and far point as 180 deg
(4) Draw several vectors v2f -- include phi_2,max]

As phi varies from 0 deg to 180 deg, v2cf spins around a semi-circle
centered at the head of vcm.  v2f will point from the tail of vcm to
the appropriate point on that semi-circle.

alpha < 1 (inelastic):

phi: 0 deg -> 180 deg
phi2: 0 deg -> phi_2,max -> 0 deg
where phi_2,max = arcsin(alpha)

alpha = 1 (elastic):  [redraw diagram]

phi: 0 deg -> 180 deg
phi2: 90 deg -> 0 deg

Note that v2f->0 as phi->0 deg and phi2->90 deg in this case.

alpha > 1 (super-elastic):  [redraw diagram]

phi: 0 deg -> 180 deg
phi2: 180 deg -> 0 deg  (all angles, including backward scattering, are possible)

Note that larger scattering angles for object 2 occur as alpha increases
(more lively collisions).  This should make sense.

Note that scattering angles for object 2 are NOT AFFECTED by changes in
beta (target/incident mass ratio).  It is true that a more massive
target will not move as quickly after the collision (|vcm| is affected),
but the scattering angle phi2 is unaffected by beta.


Note that phi2 > 90 deg (target back-scattering) requires
super-elastic collisions.

In fact, even in that case, it is difficult to conceive of target
back-scattering with simple surface contact collisions, since the
target would have to be hit "from the other side".  However,
penetrative explosive collisions (with e < 0) could result in these
angles.

This also illustrates another point: it probably makes more sense to
fix e (and not alpha) as phi (or phi') is varied, since e is much more
likely to remain constant (for given materials).  It is more difficult
to say exactly what would happen in this case, as the analysis is
quite a bit more complicated.

All of these results generalize what we found for 1-D collisions.

======================================================================

*** Elastic collisions (stationary target)
**** 1-D elastic collisions

1-D elastic collisions are a special case of the general head-on (1-D)
collisions that were discussed in earlier videos [show page].

This video is for the benefit of those who wish to review the earlier
results, as well as those who may have skipped many of the earlier
videos.

Assumptions:
  v1i, m1, m2 is given
  v2i = 0  (stationary target)
  v1f, v2f are along the same line as v1i (1-D)
  collision is elastic

Momentum conservation can be written
  m1 v1i + m2 (0) = m1 v1f + m2 v2f
==> v1i = v1f + beta v2f   (beta = m2/m1)

An elastic collision satisfies:
  |v12f| = |v12i|  (v12 = v1 - v2)
It can be shown that this condition is equivalent to
conservation of KE:
  (1/2) m1 v1i^2 = (1/2) m2 (0)^2 = (1/2) m1 v1f^2 + (1/2) m2 v2f^2
==> v1i^2 = v1f^2 + beta v2f^2

In 1-D (head-on), the condition |v12f|=|v12i| reduces to
  v1f - v2f = -(v1i - v2i) = -v1i

Putting this together with momentum conservation (2 eqns/vars) yields:
  v1f = [(1-beta)/(1+beta)] v1i = [(m1-m2)/(m1+m2)] v1i
  v2f = [2/(1+beta)] v1i = [2m1/(m1+m2)] v1i

Cases:
m1 >> m2 (very massive incident)  e.g. truck hits ball
  v1f = v1i  (incident hardly affected)
  v2f = 2 v1i
m1 > m2 (more massive incident)
  v1f = [+...] v1i  (incident continues forwards)
  v2f = [+...] v1i
m1 = m2 (equal masses)
  v1f = 0  (incident stops)
  v2f = v1i  (velocity exchange)
m1 < m2 (more massive target)
  v1f = [-...] v1i  (incident bounces backwards)
  v2f = [+...] v1i
m1 << m2 (very massive target)  ball bounces off floor/wall (earth)
  v1f = -v1i  (incident bounces back at equal speed)
  v2f = 0  (target hardly moves)

Note that in all cases, v2f is in same direction as v1i.

These conclusions are only valid if the collision is
(1) elastic
(2) 1-D (head-on)
(3) stationary target

Car example:

v1i = (+14 m/s) x-hat
m1 = 1500 kg; m2 = 2000 kg  (beta=4/3)
alpha = 1, phi = 180 deg  (head-on elastic)
vcm = v1i/(1+beta) = (+6.0m/s) x-hat
v1f = vcm (1 - alpha beta) x-hat = (-2.0 m/s) x-hat
v2f = vcm (1 + alpha) x-hat = (+12.0 m/s) x-hat

Note that v12f = v1f - v2f = (-14 m/s) x-hat = -v12i.

Kinetic energy (kJ):
            1        2      total
initial   147.0      0.0    147.0
final       3.0    144.0    147.0
change   -144.0   +144.0      0.0

======================================================================

**** 2-D elastic collisions

***** Initial setup
  m1, m2, v1i given, v2i = 0  (beta=m2/m1)
  alpha=1, but phi is arbitrary

If phi is given, we simply plug alpha=1 into the equations derived
for a general collision with a stationary target [show page]:

v1f = [v1i/(1+beta)] [(1 + alpha beta cos(phi)) x-hat + (alpha beta sin(phi)) y-hat]
    = [v1i/(1+beta)] [(1 + beta cos(phi)) x-hat + (beta sin(phi)) y-hat]
v2f = [v1i/(1+beta)] [(1 - alpha cos(phi)) x-hat - (alpha sin(phi)) y-hat]
    = [v1i/(1+beta)] [(1 - cos(phi)) x-hat - (sin(phi)) y-hat]

However, sometimes we are given phi1 or phi2 instead of phi (i.e.,
angle between v1f/v2f and v1i IN THE STATIONARY TARGET FRAME).

It might still be feasible to use the above equations in that case
(especially if phi2 is given -- with the help of trig identities one
can show that phi2 = 90deg - phi/2), but it may be more instructive to
proceed from first principles.

(a) v1f + beta v2f = v1i  [mometnum conservation]
(b) |v1f-v2f|^2 = v1i^2   [elastic condition]

Note that
|A + B|^2 = (A+B) dot (A+B)   [vectors]
          = A^2 + B^2 + 2(A dot B)
          = A^2 + B^2 + 2AB cos <(A,B)

The trick will be to express the LHS of (b) in terms of the two
velocities whose angle is given.  This can be done with the help of
(a).

======================================================================

***** Case 1: Direction of v2f is given

[Be ready with pages from "Notable features of 2-D collision -- focus on phi2"]

This case is actually quite straightforward, mainly due to the fact
that the direction of v2f is the same as the direction of J21.

To bring phi2 into the picture, we solve (a) for v1f and substitute
that result into (b):
  v1f = v1i - beta v2f
  |v1i - beta v2f - v2f|^2 = v1i^2
  |v1i - (1+beta)v2f|^2 = v1i^2
  v1i^2 + (1+beta)^2 v2f^2 - 2(1+beta) v1i v2f cos(phi2) = v1i^2
Cancel v1i^2 and move the cross term to the other side:
  (1+beta)^2 v2f^2 = 2(1+beta) v1i v2f cos(phi2)

If v2f=0, then there is no collision (pass-through elastic).

Assuming we are not interested in that solution, we divide both sides
by (1+beta)^2 v2f^1, which yields

  v2f = [2 cos(phi2)/(1+beta)] v1i
      = [2 cos(phi2) m1 / (m1+m2)] v1i

Now that we have solved for the magnitude of v2f (to go along with
the given direction), v2f is known completely.  Convert magnitude
and direction into vector components, and solve for v1f using
conservation of momentum (eqn (a)).

phi2=0 corresponds to a full head-on collision (phi=180 deg).

In fact, one can show that phi2 = 90deg - phi/2, in general.

As phi2 increases, phi decreases and the collision becomes more
"glancing".

As phi2->90 deg, v2f->0.

This is all consistent with what we said earlier, and should make
physical sense.

======================================================================

***** Ex: two cars  [draw top view now]
  m_1 = 1500 kg, m_2 = 2000 kg  (beta = 4/3)
  v_1i = 14 m/s (x-hat)
  v_2i = 0
  phi2 = 35 deg (v_2f directed 35 deg cw from +x)
  collision is known to be elastic

|v2f| = [2 cos(phi2) m1 / (m1+m2)] v1i
      = [2 cos(35deg) 1500 / (1500+2000)] (14 m/s)
      = 9.830 m/s

All velocities are in m/s...
          x        y     magnitude  direction (ccw from +x)
v_1i  +14.000    0.000    14.000       0.0
v_2f   +8.052   -5.638     9.830     -35.0
v_1f   +3.264   +7.518     8.195     +66.5+

Note: v1f = v1i - beta v2f

This same result can be obtained from the general stationary target
solution with
  alpha = 1, phi = 180 deg - 2phi2 = 110 deg

Kinetic energy (kJ):
             1        2      total
initial   147.00    0.00    147.00
final      50.37   96.63    147.00
change    -96.63  +96.63      0.0

======================================================================

***** Case 2: Direction of v1f is given

[Be ready with pages from "Notable features of 2-D collision -- focus on phi1"]

This case is more challenging.

  -v2f = (1/beta)v1f - (1/beta)v1i
  |v1f - v2f|^2 = v1i^2
  |(1+1/beta)v1f - (1/beta)v1i|^2 = v1i^2
  |(1+beta)v1f - v1i|^2 = beta^2 v1i^2
  (1+beta)^2 v1f^2 + v1i^2 - 2(1+beta) v1i v1f cos(phi1) = beta^2 v1i^2
  ... divide out v1i^2 and bring everything to one side ...
  (1+beta)^2 (v1f/v1i)^2 - 2cos(phi1)(1+beta)(v1f/v1i) + (1-beta^2) = 0
  ... factor 1-beta^2 = (1-beta)(1+beta) and cancel 1+beta
  (1+beta) (v1f/v1i)^2 - 2cos(phi1)(v1f/v1i) + (1-beta) = 0

Using the quadratic formula yields
  |v1f|/|v1i|
    = [2cos(phi1) +/- sqrt(4cos(phi1)^2 - 4(1+beta)(1-beta))] / [2(1+beta)]
    = [cos(phi1) +/- sqrt(cos(phi1)^2 - (1-beta^2))]/[1+beta]

Note that there could be zero, one, or two solutions for a given
phi1.

Negative solutions for |v1f| should generally be thrown out (a
negative solution would indicate that the actual direction of v1f is
opposite of what is given).

Once we solve for the magnitude of v1f (direction is given), v1f is
known completely.  Convert magnitude and direction into vector
components, and solve for v2f using conservation of momentum.


Case 2a: More massive target (beta > 1)

The discriminant
  D = cos(phi1)^2 - (1-beta^2) = cos(phi1)^2 + (beta^2-1)
is always positive.  Furthermore,
  |cos(phi1)| < sqrt(D)
in all cases as well.

It therefore follows that both roots are always real, with one
positive and one negative solution.  Throwing out the negative
solution yields the following:
  |v1f|/|v1i| = [cos(phi1) + sqrt(cos(phi1)^2 + beta^2 - 1)]/[1+beta]

phi1=0 represents a complete miss (phi = 0 deg).

phi1=180deg represents a full head-on collision (phi = 180 deg).

As phi1 increases from 0 to 180deg, phi also increases from 0 to
180 deg (not necessarily equal to phi1!).

This is all consistent with what we said earlier.


Case 2b: equal mass (beta = 1)

|v1f|/|v1i|
  = [cos(phi1) +/- sqrt(cos(phi1)^2 - (1-beta^2))] / [1+beta]
  = [cos(phi1) +/- |cos(phi1)|] / 2
  = 0, cos(phi1)   ...regardless of sign

Note that v1f=0 is an allowed solution, and corresponds to the full
head-on collision (phi = 180 deg).  In this case, phi1 is arbitrary,
so all possible angles phi1 admit this solution.

For the (possibly) non-zero solution...

phi1=0 represents a complete miss (phi = 0 deg).

phi1->90 deg approaches v1f=0 (the head-on collision, as a double root).

phi1 > 90deg does not lead to non-zero solutions.

Note the discontinuous behavior here -- if beta > 1, even slightly,
ALL angles to 180 deg would be allowed.  In the head-on case, if beta
is slightly larger than 1, there will be a small bounce backward.
Once beta decreases to 1 exactly, all backward angled bounces become
impossible.

Again, all consistent with our previous discussion.


Case 2c: More massive incident (beta < 1)

The discriminant
  D = cos(phi1)^2 - (1-beta^2)
will go negative, indicating no solutions, if
  |cos(phi1)| < sqrt(1-beta^2).
This may occur if
  phi1 > phi_max = arccos(sqrt(1-beta^2)) = arcsin(beta).

If phi1 > 180deg - phi_max, then cos(phi1) < 0, but D will be
positive.  However, both (real) roots will be negative (|cos(phi1)| is
greater than sqrt(D) if beta<1).

The following rules apply:
  phi1 > phi_max  => no solutions in any case
  phi1 = phi_max  => one solution (D=0; double root)
  phi1 < phi_max  => two solutions (D>0, both roots are positive)

Physically, this makes sense.  Because the incident mass is higher, it
will plow forward in all cases (there is never a solution for phi1 >=
90 deg, a generalization of the 1-D result).

There is an upper limit on the deflection angle.  That upper limit
will decrease as beta->0 --- the larger m1 is relative to m2, the
less it can be deflected by the collision.

For phi1=0, cos(phi1)=1, and there are two solutions:
  |v1f|/|v1i| = [1 +/- beta]/[1+beta]
The "+" solution represents a complete miss (phi = 0 deg); the "-"
solution represents a full head-on collision (phi = 180 deg).

For small phi1, there will similarly be two solutions: the "+"
solution (larger |v1f|) representing a glancing bounce (smaller phi);
the "-" solution (smaller |v1f|) representing a more solid, almost
head-on, collision (larger phi).

Again, this is consistent with our prior discussion: as phi
increase from 0 to 180 deg, phi1 increases from 0 to phi1_max and
then back to 0 again.

======================================================================

*** Conservation of momentum in systems
**** Momentum of a system

P_sys = sum p_i = sum m_i v_i
v_cm = sum m_i v_i / M
==> P_sys = M v_cm

The momentum of a system can be expressed entirely in terms of the
motion of its center of mass -- motion of individual particles around
the center of mass do not matter.

Split individual particle motions into motion of cm and motion around cm:
  v_i = v_cm + v_ic

Then,
  P_sys = sum m_i v_i
        = sum m_i (v_cm + v_ic)
        = sum m_i v_cm  +  sum m_i v_ic
        = M v_cm + sum m_i v_ic
        = P_cm + P_rel
  P_cm  = M v_cm = momentum due to motion of cm
  P_rel = sum m_i v_ic = momentum due to motion around cm

However,
  P_rel = sum m_i v_ic
        = sum m_i(v_i - v_cm)
        = sum m_i v_i - sum m_i v_cm
        = M v_cm - M v_cm
        = 0
Relative motion around center of mass does not contribute to overall
momentum in the system.

Momentum cannot hide.

We will be able to derive a similar partition theorem for KE:
  K_sys = K_cm + K_rel
  K_cm  = (1/2) M v_cm^2 = KE due to motion of cm
  K_rel = sum (1/2) m_i v_ic^2 = KE due to motion around cm
In this case, K_rel may not be zero.

KE can hide.

======================================================================

**** Impulse-momentum theorem for systems

Newton's second law for systems:
  F_ext = M a_cm
where F_ext = sum of all external forces acting on system (internal
forces cancel).

It follows that
  F_ext = M a_cm = M dv_cm/dt = dP_sys/dt

Integrating this over a time interval yields:
  J_ext = int(t_i->t_f) F_ext dt
        = int(t_i->t_f) dP_sys/dt
        = dP_sys (= P_sys,f - P_sys,i)

This result can also be derived directly from the impulse-momentum
theorem for individual particles, in the same manner that Newton's
Second Law for systems is derived from the same law for individual
particles:

For each particle in the system:
  dp_i = J_i,net = sum_k J_ik  [sum over all forces acting on particle i]
  J_ik = int(t_i->t_f) F_ik dt

For system of particles:
  dP_sys = sum_i dp_i = sum_i J_i,net = J_net  [sum over all forces/particles]
If we divide forces into two categories, internal and external:
  J_net = J_ext + J_int
Internal forces cancel in pairs -- so does their impulse:
  J_12 = int(t_i->t_f) F_12 dt = -int(t_i->t_f) F_21 = -J_21
and so J_int = 0

It follows that
  dP_sys = J_ext (sum of impulses due to external forces acting on system)

======================================================================

**** Conservation of momentum for isolated systems

Just as with individual particles, the impulse due to each external
force represents a transfer of momentum into the system.

The impulses associated with internal forces will transfer momentum
from one particle in the system to another, but does not represent
transfer of momentum either to or from the system overall.


Suppose there are no external forces.  Such a system is called an
isolated system.

Then F_ext = 0, and so J_ext = 0 during any time interval.  It follows
that dP_sys = 0 and so P_sys,f = P_sys,i.  Momentum is conserved for
such systems.

This can equally well be derived this way:
  dP_sys/dt = F_ext = 0 ==> P_sys = constant

Momentum can shift around within the system as a result of internal
forces, but the overall momentum of the system cannot change.

Note that J12 = -J21 for internal forces (Newton's 3rd Law), and so
any changes in momentum of particle 1 in the system will be offset
by changes in momentum of particle 2 in the system.


When two objects collide (assuming that non-collision forces can be
neglected during the collision), the two-object system can be
approximated as isolated during the collision time interval, and thus
the total momentum of the two objects does not change during the
collision (i.e., momentum is conserved during the collision).

This is normally how conservation of momentum for two-body collisions
is proved.

Note that the internal "collision" force can result in significant
momentum transfer BETWEEN the two objects, even as total momentum
does not change.

The result we just derived extends this conservation principle to
multi-particle systems.


Note that the only stipulation on momentum conservation is the
isolated system (no interactions with the outside world).  Momentum
conservation is a universal conservation principle.

When one gets more sophisticated in addressing "action at a distance"
type forces (e.g., gravity, electromagnetic, etc), one will find that
Newton's 3rd Law (as we have currently stated it) will fail.  However,
the fields that mediate such interactions can, themselves, carry
momentum, and one finds that the momentum conservation principle
continues to be upheld.

======================================================================
NOT-YET

**** Multiparticle explosion example [come up with one]

======================================================================

======================================================================

**** Bag of sand dropped on cart example

[draw diagram; indicate initial velocities]
  m_1 = 20 kg  (cart)
  m_2 = 7.0 kg  (sand)
  v_1i = (+15 m/s) x-hat
  v_2i = 0 x-hat - (5.6 m/s) y-hat   (dropped from 1.6 m)
  v_f = v_1f = v_2f = ?

Naive solution:

m_1 v_1i + m_2 v_2i = m_1 v_1f + m_2 v_2f = (m_1 + m_2) v_f
==> v_f = [(300 kg m/s)x-hat - (39.2 kg m/s)y-hat] / (27 kg)
        = (11.11 x-hat - 1.46 y-hat) m/s

The problem?  v_f cannot have a negative y-component, because the
cart cannot pass through the floor.  In fact, the cart (and the bag
of sand) must both travel together in the x-direction.

This is an example of a collision between two objects for which
external forces are NOT negligible during the collision.  We therefore
cannot assume that momentum is conserved.

How do we handle this?  We can still apply the impulse-momentum theorem
for systems:
  P_sys,f - P_sys,i = J_ext

Force diagrams:
Cart(1):
EXT  m1 g (down)
     N_1G (up)
INT  N_12 (down)
     f_k12 (left)
Sand(2):
EXT  m2 g (down)
INT  N_21 (up)
     f_k21 (right)

The issue is with N_1G, an external force.  Because the bag of sand is
pushing down very strongly on the cart (N_12), the ground must push up
with the same order of magnitude force (N_1G).  N_1G is not
negligible, and thus J_N1G cannot be assumed to be zero.

Note, however, that all external forces are vertical (up/down).  It
follows that J_ext,x is zero, even if J_ext,y is not zero.

Thus
  P_sys,fx - P_sys,ix = J_ext,x = 0  ==> P_sys,fx = P_sys,ix
  P_sys,fy - P_sys,iy = J_ext,y != 0  ==> no such conservation

Applying momentum conservation in the x-direction yields
  v_fx = [m_1 v_1ix + m_2 v_2ix] / [m_1 + m_2] = 11.11 m/s
We don't have momentum conservation in the y-direction, but we don't
need it -- v_fy=0 because the motion of the cart + bag is constrained
that way.
    v_fy = 0  (constraint)

==> v_f = (11.11 m/s) x-hat

Note that f_k12 and f_k21 (internal collision forces) act in the
x-direction, and act in a way to reduce the relative velocity between
the cart and the sand to zero.  f_k12 acts to the left, and slows the
cart.  f_k21 acts to the right, and causes the bag to move with the
cart.

N_12 and N_21 (also internal collision forces) act in the y-direction,
and act in a way that stops the bag of sand from falling.  Without the
floor in the way, the cart would be pushed downwards, and the two would
move off with a final velocity of (-1.46 m/s) in the y-direction.

However, the ground IS in the way, and N_1G acts as well.  We can
calculate J_N1G using the impulse-momentum theorem:
  J_N1G,y = J_ext,y = P_sys,fy - P_sys,iy = (0 kg m/s) - (-39.2 kg m/s)
          = +39.2 N s

Note that impulse due to gravity (on the cart and the bag) should
technically be included in J_ext, but they are negligible (since dt
is short), and can be safely neglected.

======================================================================
